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ABSTRACT 

We study the "normal" decay phase of the X-ray afterglows of gamma-ray 
bursts (GRBs), which follows the shallow decay phase, using the events simulta- 
neously observed in the R-band. The classical external shock model — in which 
neither the delayed energy injection nor time-dependency of shock micro-physics 
is considered — shows that the decay indices of the X-ray and R-band light 
curves, ax and ao, obey a certain relation, and that in particular, ao — ax 
should be larger than —1/4 unless the ambient density increases with the dis- 
tance from the central engine. For our selected 14 samples, we have found that 
4 events violate the limit at more than the 3a level, so that a fraction of events 
are outliers of the classical external shock model at the "normal" decay phase. 
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Introduction 



Gamma-ray bursts (GRBs) consist of two phases: prompt GRB emission and subsequent 
afterglows. How long the prompt GRB emission lasts and when the transition from the 
prompt GRB to the afterglow occurs have been long-standing problems. These problems 
are tightly related to the mechanism of the central engine of GRBs. The Swift satellite has 
brought us early, dense, and detailed data on the afterglows of GRBs in various observation 
bands. Now, we are entering the era of multi- wavelength observations; especially optical and 
X-ray bands, which tell us some hints for answering the problems. 

Contrary to the expectation in the pie- Swift era, ^wz/i^ X- Ray Telescope (XRT) data have 
revea l ed complex temporal behavior of the X-ray afterglow (Burrows et ahlboOS : Tagliaferri et al 



200,4 iNousek et al.l l2006l : lO'Brien et al.ll2006al : IWillingale et al.l 120071 ). Initially, it decays 



yery steeply, whose inost popular interpretat i on is the tail emissio n of the prompt GRB 



( Kumar &: Panaitescu 200ol : Zhang et al 
bilities have been proposed (e.g. 



2006 



Yamazaki et al.ll2006l ). although other possi- 



Zhang et al.ll2007l ) . At several hundreds of seconds after the 



Toma et al. 


2006; 


loka et al. 


2006; 


Zhang 


2007) 



(e.g., 

the X-rays subsequently decays with the decay index usually steeper than unity, which was 
expected in the wg- S wift era. This d ecay behavior can be well explained by the classical 
external shock model ( ISari et al.l Il998l ). in which neither the delayed energy injection nor 
time-dependency of shock micro-physics is considered. Hence this phase is sometimes called 
the normal decay phase. 

However, as the number of the X-ray observations increases, it is getting suspicious 
that the normal decay phase arises from the external shock. In the steep and the shallow 
decay phases, the X-ray light curves someti mes possess large bumps, called the X-ray flares 
( IChincarini et al.l 120071: iFalcone et al.l 120071) . and/or dips that cannot be explained by the 
external shock model (lloka et al.ll2005l ). Furthermore, for an extreme example, GRB 070110 
showed a rather complex X-ray afterglow with a s udden drop at ~ 2 x 10^ sec after the 
burst trigger as the end of the shallow decay phase (jTroja et al.l 120071 ). These observational 
facts may tell us that the steep and the shallow decay phases are likely due to late internal 
dissipation of the energy produced by the long-acting central engine. On the ot her hand, the 



X-ray spectrum remains unchanged across the shallow-to-normal transition (INousek et al. 



20061 ). which may imply that the shallow and the normal decay phases are of the same origin. 
Therefore, it might be that the normal phase comes from the internal energy dissipation. 
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The observed optical afterglow is also com plicated and in an early epoch (< 10^ sec), 
there is a diversity (jZhangll2007l : iDoi et al.ll2007l ). On the other hand, it was found in the pre- 
Swift era that for almost all events, the behavior at > 0.1 day after the burst could be well 



expla ined by the c/asszca/ external shock model (e.g.. |Panaitescu fc Kumarll200ll : IUrata et al. 



20031). although sorne events showed c omplex light curve s with dips and/or bumps (e.g. 



Holland et al. 



20031 : iLipkin et al.l 12004 lUrata et al.ll2007bl ). This epoch corresponds to the 



normal decay phase of the X-ray afterglow. Those previous studies are mainly based on the 
optical bands, because the X-ray observation was sparse at that time. In the Swift era, we 
are starting to have the simultaneous optical and X-ray afterglow data in the epoch > 0.1 
day after the burst thanks to the rapid and the dense X-ray observation by the XRT. 

In this Letter, we study the normal decay phase of the X-ray afterglows simultaneously 
observed in optical R-bands, and investigate whether it is consistent with the classical ex- 
ternal shock model or not. We perform a simple test using the optical and the X-ray decay 
indices, ao and ax, where we use a notation, F^, oc t~"z/~^. For example, in the classical 
external shock model with uniform ISM, they are related to the power-law index of the 
electron distribution, p{> 2), as ao = 3(p — l)/4 and ax = (3p — 2)/4, r espectively, sinc e 
the cooling frequency Uc usually lies between the optical and X-ray bands (ISari et al.lll998l ). 
Eliminating p, w e obtain an — ax = — 1/4. Similarly, for the wind environment, we derive 
ao — ttx = 1/4 (IChevarier fc Lil 120001). Thes e relations between ao and ax are also valid 
in the case of 1 < p < 2 ( jPai fc Cheng) l200ll ). Therefore, through the relation between ao 
and ax, one can test the classical ext ernal shock model. In the pre-Swift era, similar study 
has been done for BeppoSAX GRBs (jPe Pasquale et al.ll2006bl ). However, compared with 
the Swift GRBs, their X-ray data were not well enough to identify the normal decay phase 
and to determine the decay index with small uncertainties. We can now obtain more dense 
X-ray and optical data and can determine ao and ax with much less ambiguity. Finally 
it is noted that in this Letter, we do not consider the spectral indices, /Sq and /3x, because 
they have at present large uncertainties; /?o fairly depends on the assumed dust model, and 
the low X-ray flux at the epoch we are interested in makes us difficult to constrain /3x with 
precision which we need to test the model. 



2. Decay indices of X-ray and R-band afterglow in the normal decay phase 

We consider long GRBs that are followed-up by Swift XRT from the beginning of 2005 
to the end of 2006. The Swift XRT data are systematically analyzed using our pipeline script. 
The cleaned event data of the Window Timing (WT) and the Photon Counting (PC) mode 
from the Swift Science Data Center (SDC) are used in the whole process. Although both WT 
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and PC mode data are processed in the pipeline, hereafter, we are only focusing on the process 
of the PC mode data. The search of the X-ray afterglow counterpart, a construction of the 
X-ray light curve, and a fitting process of the X-ray light curve and spectra are performed 
automatically using the standard XRT softwares and calibration database (HEASoft 6.2 and 
CALDB 20070531). The source region is selected as a circle of 47" radius. The background 
region is an annulus of an outer radius of 150" and an inner radius of 70" excluding the 
background X-ray sources detected by ximage in circle region of 47" radius. The light curve 
is binned based on the number of photons required to meet at least 5(j (Sakamoto et al. 
2007). We select the samples of the X-ray afterglows which have a smooth transition from 
the shallow to the normal decay phases at > 10'^ sec. Samples with X-ray flares have been 
excluded. Then, we find the start time of the normal decay phase of the X-ray afterglow 
(ax ^ 1); and extract events in which well-sampled R-band light curves are available during 
the normal decay phase. 

The light curve data in the R-band ar e published in lite ratures or observed by the East 
Asian ORB Follow-up Network {EAFON; lUrata et all l2005[ ll^ and the KANATA telescope. 
The R-band data taken by us are processed as in the following. A standard routine, including 
bias subtraction, dark subtraction, and fiat-fielding corrections with appropriate calibration 
data is employed to process the data using IRAF. Flux calibrations are performed usin g the 
APPHOT package in IRAF, referring to the standard stars suggested by lLandoltl (Il992l ). For 
each data set, the one-dimensional aperture size is set to 4 times as large as the full- width 
at half maximum of the objects. The magnitude of error for each optical image is estimated 
as cr^ = cTpjj + CTgyg, where Uph represents the photometric errors for each afterglow, estimated 
from the output of IRAF PHOT, and (Xsys is the photometric calibration error estimated 
by comparing our instrumental magnitudes. When we combine data which are obtained 



Huang et al. 


2007; 


Urata et al. 


2007a. 


2003) 



and yield realistic light curves. 

There are 14 GRBs which have a good coverage with both X-ray and optical bands at 
the normal decay phase. Among them, optical data of 11 events have been already published 
in literatures. For unpublished data obtained by EAFON, detailed light curves in the X-ray 
and the optical bands are presented elsewhere (Urata et al. 2007, in preparation). For those 
samples, we identify the normal decay phase that is well described by a single power-law 
decay model and derive ax- During the phase, we find that the optical light curves are well 
fitted with a single power-law model in the time interval shown in Tabled] in which the decay 



^In this paper, the sa mples are mainly ta ken using Kiso 1.05m Schmit ed telescope j Urata et al. 20031) , 
Lulin One-meter telescoDe (|Huang et al.ll2005| ) and Xinglong 0.8m telescoDe ljDeng et al.ll2007l) . 
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index, ao, is determined. All results are summarized in Tabled! Figure [T] shows ao as a 
function of ax, while Figure [2] shows the value of ao — ax for each event. The quoted errors 
in this Letter are at the 1 a confidence level. 



Results and Discussion 



Let us consider the case of the minimum frequency smaller than the R-band frequency 
(^m < ^r), which is a reasonable assumption for several bright bursts in the pre- Swift era. 
In this case, the spectral index of the optical afterglow is positive, /3n > 0, which is consistent 
with the previous observational results (see Table 2 of lKann et al.ll2006l ). The decay and the 
spectral indices are calculated as shown in Table [2]by the classical external shock model with 
ambient matter density dependent on the radius, n cx r~* where we assume s > 0. Since 
the Lorentz factor of the relativistically expanding shell evolves with the observer time as 



F oc t" 
derive 



s < 3 is needed in order for the shell to decelerate. If i^m < i^r < i^c < i^x, we 



1 



-2s 



which is valid for 1 < p < 2 or 2 < p, so that ao — ax ranges between —1/4 and 5/4 if 
< s < 3. For the cases of i^m < J^r < t^x < t'c or i^m < < J^r < t^x, «o — «x should be zero. 
As can be seen in Fig. |2l among 14 events considered in this paper, 4 events (GRB 050319, 
050401, 060206, 060323) are below the line ao — c^x = —1/4 at more than the 3a level, 
so that a fraction of bursts are outliers of the classical external shock model at the normal 



decay phase. iDe Pasquale et al.l (l2006bl ) have performed similar study and found two out 
of 12 events have ao — ax significantly below —1/4, which was roughly consistent with our 
result (see Table 5 of their paper) . 



Liang et al.l (120071 ) studied the ax^Px relation of the normal decay phase and found 
that there are several outliers of the classical external shock model. Their outliers have 
large ax > 2. In our sample, however, outliers of ao^ax relation exists even if their ax 
of around 1.5, and their ax~/3x relat ions are consistent with the classical external shock 
model (see Fig. 5 of iLiang et al.l 120071 ). This fact, therefore, strengthens the importance of 
the multi-wavelength studies at the normal decay phase to test the classical external shock 
model. 

There are several pos sibilities leading to ao — ax < — 1/4. One is to consider s < 



case (e.g. lYost et al.ll2003l ). If s < —4, then ao — ax ^ —0.5, however, there is no theo- 
retical reason to consider such a steeply rising profile. Another is to consider the delayed 
energy injection (IRees fc Meszaroslll998l ) and/or time- variable shock micro-physics parame- 
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ters (lYost et al.ll2003l ). Here, we conside r Um < j^r < Vc. <yx for simplicity. The generalized 
forms of ao and ax are then derived by iPanaitescu et al.l (120061 ) based on the assumptions 
E{> r) cx r~^, Eb o(. r~^, oc r~*, and n oc r~*, where s < 3. Then, from their derived 
formula, we obtain 



1 
4 



8 -2s 
3-s 



4(e + 8 - 2s) V 4 - s 



4 -3s 



e + 3b 



which is independent of i and p. We find that ao — «x < ~l/4 is achieved if e + 36 > for 



the uniform ISM case (s = 0), or if 6 — e > 8/3 for the win d medium case 



s = 2). Note that 



these cases have been discussed for the pre-Swift GRBs ( jPiro et al.l 1 19981 : lYost et al.l 12003 
Corsi et al.l 120051 ). 



Another possibility to explain the outliers may be that the X-ray flares superposing 
on the X-ray afterglow could steep the apparent decay index of the X-ray. X-ray flares are 
usually more active in the initial phase, so that they may enhance the early X-ray flux. In 
this case the X-ray flare should not be spiky but relatively smooth, and the late afterglow is 
just an ordinary afterglow. 

Although the external shock model is still viable, the afterglow emissions of outliers 
may be capabl e of th e internal shock origin. Such a possibility has been proposed by 
Ghisellini et al.l (120071 ). Then, the optical and X-ray emission in the late phase are of differ- 
ent origins. It is also possible in this model that a chromatic break occurs at ~ 1 day after 
the burst, which was beli eved to be achromatic in the we-Swift era and to be cau sed by 
the jet coUimation effects (IPanaitescu et al.ll2006l: iHuang et al.ll2007l : ISato et al.ll2007r ). Or a 
cannonball model may account for our outliers (jPado et al.ll2007l ). 
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CYopt 




Fig. 1. — The R-band decay index ao as a function of the X-ray decay index ax in the 
normal decay phase. The classical external shock model predicts uq — ay^ = —1/4 (solid 
line) and 1/4 (dashed line) for the uniform ISM (s = 0) and for the wind medium (s = 2) 
cases, respectively. 
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Table 1. X-ray and optical temporal decay indices during X-ray normal decay phase. 
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.3 


X 


10* 


- 2, 


.5 


X 


10* 


1, 


..39 


± 


0. 


,08 


0, 


.86 




0, 


.03 


-0.53 ±0.09 


(1), 


(11), (12) 


060206.2 


2, 


.3 


X 


10* 


- 5, 


.4 


X 


10^ 


2, 


.5 


X 


10* 


- 2, 


.0 


X 


10" 


1, 


.39 


± 


0. 


.08 


1. 


.41 




0. 


.02 


0.02 ± 0.08 


(1), 


(11), (12) 


060323 


1. 


.1 


X 


103 


- 2. 


.1 


X 


10^ 


1, 


.2 


X 


103 


- 3, 


.0 


X 


10^^ 


1, 


.38 


± 


0. 


.10 


0. 


.70 


± 


0. 


.05 


-0.68 ±0.11 


(1) 




060526 


1. 


.8 


X 


10* 


- 4. 


,2 


X 


10^ 


2, 


.0 


X 


10* 


- 3, 


.2 


X 


10* 


1, 


.60 


± 


0. 


.13 


1. 


.17 


± 


0. 


.07 


-0.43 ±0.15 


(13) 




060605 


5. 


,2 


X 


103 


- 2. 


,7 


X 


10* 


2, 


.0 


X 


10* 


- 2, 


.3 


X 


10* 


1, 


.45 


± 


0. 


.06 


1, 


.41 


± 


0, 


.44 


-0.04 ± 0.44 


(1) 




061121.1 


2. 


.1 


X 


103 


- 1, 


,7 


X 


10* 


4, 


.7 


X 


103 


- 1, 


.5 


X 


10* 


0, 


.98 


± 


0. 


,04 


0, 


.98 


± 


0, 


.05 


0.00 ± 0.06 


(14) 




061121.2 


1. 


.7 


X 


10* 


- 3, 


,5 


X 


10^ 


7, 


.2 


X 


10* 


- 3, 


.3 


X 


10® 


1, 


.46 


± 


0. 


.03 


1, 


.54 




0, 


.07 


0.08 ± 0.08 


(14) 


, (15) 



^Thc nomal decay phase is identified in the Swift XRT data. The value of ax is determined in this period. Time zero is taJcen 
as the burst trigger time. 

''The period when the optical data was talcen during the normal decay phase. The value of ao is determined in this epoch. 

■^References for optical data. (1) EAFON (for specific individual events, e.g. Huang et al. (2007), Deng et al. (2007)); (2) De 
Pasquale et al. (2006b); (3) de Ugarte Postigo et al. (2007); (4) Klotz et al. (2005); (5) Delia Valle et al. (2006); (6) Antonelli 
et al. (2006); (7) Rykoff et al. (2006); (8) Cenko et al. (2006); (9) SoUerman et al. (2007); (10) Yost et al. (2007); (11) Wozniak 
et al. (2006); (12) Stanek et al. (2007); (13) Dai et al. (2007); (14) Uemura et al. (2007); (15) Halpern et al. (2007) 
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Table 2. Spectral and decay indices {F,y oct °'v ^) predicted by the classical external 

shock model. 



l<p<2 2<p 

a (5 a (5 



s-2 _ 1 
4-s 3 

3p— 1 _|_ s p— 1 
4 8-2s 2 
3p-2 p 
4 2 



m < i/ < f c 



(4s-3)p-2(s+3) 
3(p-l)(8-2s) 
(3-s)p+(6+s) 

2(8-2s) 
(3-5)p+2(5-s) 
2(8-2s) 



_ 1 
3 

p-1 
2 

2 



^Note: For the case of the spherical expansion, 
slow cooling, and the ambient density profile given by 
n oc r~*, where < s < 3. The break frequency 
i^m evolves with time as oc t~^/'^ for 2 < p while 

(3 — s)p+6 — 5 

f m oc t (p-i)(8-2s) for 1 < p < 2, and scales as 

3a-4 

OC ^8-2= regardless of p. 



